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ABSTRACT 

This  study  proposes  to  analyze  the  correlation  sensitivity  of 
si>atial  Wiener  spectra  derived  from  polarized  (Stokes'  pareuMter) 
imagery  of  both  synthetic  and  measured  waterwave  scenes.  The  intent 
of  this  study  is  to  evaluate  the  utility  of  polarized  imagery  in 
the  extraction  of  spatial  information  content  (to  include  texture) 
from  scenes  containing  waterwaves,  one  of  many  spatially  correlated 
natural  features  that  can  be  described  by  Wiener  spectra. 

This  study  will  first  enhance  an  existing  simulation  model  [North, 
1989]  to  support  the  comparison  of  the  spatial  performance  of  an 
imaging  system  through  the  creation  of  synthetic  polarized  images 
of  waterwave  scenes  and  their  spatial  spectra.  Second,  this  model 
will  be  tailored  to  a  specific  imaging  system:  differentially 
polarized  4-lens  35-mm  photography.  Third,  a  parametric  exploration 
of  the  model  will  be  used  to  predict  the  polarized  imaging 
conditions  that  optimize  spatial  information  content.  Fourth,  a 
polarized  4-lens  camera  system  will  collect  representative  imagery 
containing  waterwave  scenes  under  natural  conditions.  Finally,  the 
synthetic  and  measured  polarized  image  spectra  will  be  analyzed  euid 
compared  in  order  to  validate  the  simulation  model,  calibrate  the 
imaging  system,  and  ultimately  provide  the  basis  for  the  utility 
evaluation. 

The  results  of  this  study  should  provide  a  refined  understanding  of 
the  contribution  that  sensed  polarization  can  provide  in  the 
analysis  of  spatial  information  content  from  imagery,  through  the 
detailed  examination  of  its  effect  on  natural  waterwave  surfaces. 
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1.0 


PURPOSB 


The  purpose  of  this  study  is  to  evaluate  the  utility  of  polarized 
(Stokes'  paraseter)  imagery  in  the  extraction  of  spatial 
information  content  from  scenes  of  waterwaves,  one  of  many 
spatially  correlated  natural  features  that  can  be  described  by 
Wiener  spectra.  The  primary  method  of  evaluation  is  correlation 
sensitivity  analysis  of  the  synoptic  spatial  spectra  derived  from 
the  polarized  imagery. 

This  investigation  will  be  executed  in  three  parts: 

1)  synthesize  a  reasonable  model  which  can  consider  the  effects  of 
resolving  differentially  polarized  spatial  image  spectra  of 
waterwave  scenes  under  various  illumination,  scene,  and  sensor 
geometries;  tailor  the  model  for  a  specific  imaging  system  - 
differentially  polarized  4-lens  35-mm  photography;  and  exercise  the 
model  to  define  the  parametric  dependencies  for  correlation  between 
the  actual  (modeled)  waterwave  slope  spectra  and  their 
reconstructions  from  the  polarized  radiance  distributions  captured 
by  imagery, 

2)  execute  a  controlled  experiment:  photograph  representative 
waterwave  surfaces  under  natural  illiuination  conditions  by  a 
differentially  polarized  4-lens  35-mm  camera  system,  using  the 
model  predictions  to  specify  the  range  of  imaging  conditions  that 
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should  enhance  the  correlation  sensitivity  analysis,  and 

3)  analyze  and  conpare  both  the  synthetic  and  measured  polarized 
image  spectra;  attempt  to  validate  the  simulation  model;  attempt  to 
calibrate  the  imaging  system;  and  provide  a  utility  evaluation. 

The  prime  motivation  for  this  study  is  to  extend  an  existing 
simulation  model  of  one  spatially  correlated  natural  feature 
(waterwaves)  as  a  basis  for  the  study  of  other  natural  background 
and  clutter  features  (definable  by  Wiener  spectra)  under  polarized 
imaging  conditions.  However,  this  simulation  model  is  unvalidated, 
hence  a  second  motivation  is  to  calibrate  the  model  through  the 
analysis  of,  and  comparison  with,  actual  polarized  image  data.  A 
third  motivation  is  to  consider  the  potential  utility  of  polarized 
waterwave  scenes  as  natural  opportunistic  targets  for  synoptic 
image  calibration.  This  motivation  is  based  on  several  unique 
aspects  of  water  and  waterwaves  that  are  enumerated  below.  A 
fourth,  and  final,  motivation  is  to  demonstrate  the  feasibility  of 
extending  this  analysis  to  more  complex  imaging  systems,  e.g., 
synthetic  aperture  radar  (SAR)  polarimetry,  using  simpler  optical 
techniques  for  this  study. 

There  are  several  reasons  for  selecting  natural  waterwave 
surfaces  as  the  object  of  polarized  image  spectzxun  analysis: 

1)  Terrestrial  surface  water  is  ubiquitous;  fully  three 
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quarters  of  the  Earth's  surface  is  covered  with  water.  For  airborne 
and  spaceborne  sensors,  this  represents  a  large  window  of 
opportunity  to  image  this  target. 

2)  In  comparison  with  other  terrestrial  surface  features, 
water  surfaces  are  optically  homogeneous;  only  two  radiance 
functions,  one  for  surface  reflection  and  one  for  upwelling 
subsurface  refraction,  are  required  to  describe  the  unique  radiance 
contribution  of  the  imaged  water  scene.  However,  in  the  red  (and  at 
longer  wavelengths ) ,  water  is  increasingly  opaque  to  the  point  that 
only  surface  reflection  (and/or  emission  in  the  infrared)  requires 
description. 

3)  Except  for  the  rare  instances  of  sustained  windlessness 
idiere  water  is  smooth  and  specular,  its  surface  is  disturbed  by  the 
addition  of  directional  friction  energy  applied  by  wind  passing 
over  it.  Up  to  a  defined  maximum  wind  velocity,  the  disturbed 
surface  remains  analytic  (infinitely  differentiable)  and  can  be 
described  as  a  quasi-stationary,  pseudo~Gaussian  process;  the 
resulting  power  spectrum  approximations  for  elevation,  slope,  and 
cizrvature  can  be  analyzed  by  Fourier  methods  [Kinsman,  1965 ] . 

4)  As  an  extension  of  3)  above,  waterwaves  can  be  imaged  over 
several  orders  of  spatial  frequency.  Also,  the  well-behaved  surface 
statistics  of  waterwave  surfaces  (e.g.,  mean  slope  =  0*,  3-sigma 
slope  «  30*)  allow  for  a  broad  range  of  off-vertical  imaging 
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geoaetries  without  having  to  consider  the  complex  effects  of 
surface  obsciiration. 

5)  Also  as  an  extension  of  3),  the  assumption  of  quasi-spatial 
and  quasi-temporal  invariance  of  the  surface  spectra  allows  for 
considerable  tolerance  in  both  image  collection  and  processing: 
precise  pointing  to  an  exact  ground  reference  point  by  multiple 
sensors  is  not  critical;  precise  collection  timing  by  multiple 
sensors  is  not  critical;  and  image  registration  is  not  required. 

6)  Finally,  the  natural  air-water  interface  represents  a 
strong  polarizing  dielectric.  High  contrasts  within  emd  between 
polarized  images  should  be  readily  detectable  in  a  significant 
range  around  mean-Br«Nrster-angle  imaging  geometries. 

In  sum,  this  study  will  exploit  one  of  the  more  analytically 
straightforward  terrestrial  surface  features  for  the  evaluation  of 
polarized  (Stokes'  parameter)  imagery  in  the  extraction  of  texture 
information  from  scenes  of  spatially  correlated  natural  features 
that  can  be  described  by  Wiener  spectra. 
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2.0 


LITBRATDRB  REVIEW 


2.1  laoig*  Analysis  of  Watervave  Surfaces 

In  1925,  Schumacher  made  near-horizontal  stereo-pairs  from  a 
ship  with  the  intent  of  measuring  the  variability  of  wave  heights. 
The  utility  of  this  method  was  severely  limited  due  to  many 
factors:  1)  the  camera  baseline  was  restricted  to  the  length  of  the 
ship;  2)  waves  in  the  foreground  obstructed  %iaves  in  the 
background;  3)  backsides  of  waves  were  not  visible;  and  4)  there 
was  a  lack  of  'ground'  control  on  the  open  seas  for  height 
determination  -  the  errors  are  especially  pronounced  from  an 
oblique  perspective  [Pos,1988].  This  experiment  represents  a 
limiting  case  for  watermve  surface  measurement  via  image  analysis. 

In  1933,  Hulburt  [1934]  made  polarized  and  unpolarized  oblique 
photographs  of  sun  glitter  on  sea  waves  to  measure  the  polarization 
of  light  at  sea  with  respect  to  surface  roughness,  sun  angle,  and 
weather  conditions.  Because  the  widths  of  glitter  patterns 
correlate  to  the  maximum  slope  of  the  sea  surface,  Hulburt  was  able 
to  demonstrate  that  waves  in  the  North  Atlantic  varied  from  15* 
inclination  when  winds  were  blowing  at  3  knots  up  to  25* 
inclination  at  18  Jcnots. 

m 

Sawyer  [1949]  mounted  a  Sonne  strip  camera  on  a  fast  low- 
flying  airplane  to  photograph  narrow  strips  of  the  sea  surface  and 
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neasure  the  directional  spectrun  of  the  waves.  The  field  of  view 
was  too  narrow  to  capture  significant  anounts  of  surface  data 
orthogonal  to  the  flightline.  The  accuracy  of  the  spectral  estinate 
decreased  with  the  angle  from  the  flightline. 

Also  in  1949,  Barber  [1949,1954]  analyzed  single  photographs 
of  sea  surfaces  to  detemine  wave  direction,  but  he  was  unable  to 
determine  the  two-dimensional  spectrum  due  to  the  computing 
limitations  of  the  time.  At  this  point,  it  became  apparent  that  an 
essential  requirement  for  future  analysis  of  waterwave  surfaces  was 
to  have  near-nadir,  high  resolution,  high-contrast  images  covering 
large  areas.  The  intent  was  to  photographically  capture  significant 
information  about  the  largest  range  of  spatial  frequency  components 
without  perspective  distortion  or  hidden  surface  detail.  It  also 
became  apparent  that  the  analysis  of  large  amounts  of  spatial 
information  required  the  power  of  a  digital  computer. 

In  1951,  cox  and  Munk  [  1954a, b]  measured  the  wave-slope 
distributions  of  the  sea  surface  from  aerial  photographs  of  sun 
glitter  patterns.  They  computed  the  distribution  from  the  measured 
variation  of  radiance  within  a  glitter  pattern  instead  of  computing 
maxima  from  the  pattern  boundaries  as  done  by  Hulburt.  Four  cameras 
were  flown  from  a  single  airplane  at  altitudes  of  2000  feet,  with 
two  used  as  imagers  and  two  used  as  radiometers.  Their  image 
analysis  was  quite  sophisticated;  it  accounted  for  sun  diameter, 
angular  reflectivity,  lens  falloff  (vignetting),  film  sensitivity. 
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exposure  calibration  (sensitoaetry) ,  and  ultiaately  provided  a 
first-order  relationship  between  filn  density  and  directional  wave- 
slope  probability. 

In  1953,  Schooley  [1954]  performed  a  sinplified  version  of  the 
Cox  &  MunX  experiment  by  taking  flash  photographs  of  a  river 
surface  from  a  45-foot  bridge  elevation  at  night.  The  main 
limitation  of  his  experiment  %#as  the  probable  inhomogeneity  of  the 
water  surface  due  to  limited  fetch  and  the  presence  of  wave- 
refracting  obstacles  in  the  water. 

In  1954,  medium-altitude  (3000  feet)  stereophotography  was 
employed  by  Marks  and  Ronne  [1955]  to  generate  stereopairs  of  sea 
surfaces.  Two  airplanes  carried  radio-synchronized  cameras  and  a 
surface  ship  acted  as  'ground'  control  in  the  photographs. 
Elevations  were  photogrammetrically  measured  at  discrete  points  and 
the  sampled  elevation  array  was  then  autocorrelated  (the  sampling 
distance  determined  the  desired  spatial  resolution) .  This 
experiment  marks  the  first  recorded  use  of  a  digital  computer  to 
calculate  the  directional  2D  spectra  of  waterwaves.  The  work  of 
Cote  et  al.  [1960]  enhanced  this  basic  technique.  More  recent 
stereophotogrammetric  efforts  include  Holthuijsen  [  1983a, b]  and  Pos 
et  al.  [1988].  Elements  of  this  later  work  include  methods  to 
render  the  water  opaque  so  that  a  more  exact  calculation  of  the 
height  field  can  be  made. 
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During  the  1950s  and  early  1960s,  Longuet-Higgins  [1952-1962; 
Cartwright  &  Longuet-Higggins,1956]  elaborated  on  the  results  of 
Cox  and  Munk  to  fomulate  the  statistical  theory  of  patterns, 
paths,  niunber,  frequency,  and  distributions  of  specular  reflection 
points  on  randomly  moving  surfaces.  Stilwell  [1969;  Stilwell  & 
Pilon,1974]  correlated  the  statistics  of  sea-surface  images  to  the 
wave-slope  statistics  of  the  actual  sea  surface.  Under  the 
assumptions  of  uniform  sky  radiance,  optimized  viewing  geometry, 
and  small  surface  slopes,  Stilwell  derived  the  relationship  between 
the  film  transmittance  of  an  imaged  surface  point  and  the  range 
component  of  the  wave  slope  at  that  surface  point.  He  further 
demonstrated  a  linear! zable  relationship  between  the  spatial  image 
spectrum  and  the  imaged  surface  slope  spectrum.  Kasevitch  [1975] 
extended  Stilwell 's  model  to  second  order  to  develop  an 
optimization  criterion  for  the  relationship.  Chapman  and  Irani 
[1981]  took  this  work  one  step  further  by  creating  a  synthetic 
model  and  executing  a  limited  quantification  of  the  error  magnitues 
associated  with  the  parametric  dependence  of  this  linear  modc^l. 

Sheres  [1980]  developed  a  novel  technique  for  remotely  sensing 
surface-flow  velocities  based  on  imagery  of  monochromatic 
wavetrains  of  known  frequency  (such  as  those  generated  by  a  motor 
boat)  propagating  over  the  region  of  interest.  His  work 
demonstrated  that  the  wavelength  and  direction  of  two  different 
wavetrains  produced  all  the  information  required  to  calculate 
surface  flows.  Gotwols  and  Irani  [1980]  developed  a  similar 
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technique  to  determine  the  phase  velocity  of  short  gravity  waves. 


Exotic  sensors  using  LASER  [Palm  et  ai.,1977;  Schau,1978; 
Abshire  &  McGarry,1987]  and  LIDAR  [Weinman, 1988]  have  been  used  to 
extract  directional  spectra  and  surface  backscatter  data  at  higher 
wavenumbers  (i.e.,  the  capillary  wave  regime).  Synthetic  Aperture 
Radar  (SAR)  iamgery  has  been  used  to  estimate  spectra,  phase 
velocities,  and  propagation  directions  at  lower  wavenumbers  (i.e., 
the  gravity  wave  regime)  [Monaldo  &  Lyze'iga,1986;  Nonaldo  & 
Kasevich,1982;  Carlson, 1984 ] .  Also,  Long  Wave  Infrared  (LWIR) 
sensors  have  been  used  to  calculate  spatial  spectra  of  ocean- 
surface  temperature  [Saunders, 1967, 1968;  McLeish,1970] . 

Lybanon  [1985]  reported  on  the  implementation  of  an  automated 
image-analysis  system  by  the  U.S.  Naval  CX:ean  Research  & 
Development  Activity  (NORDA) .  The  Interactive  Digital  Satellite 
Image  Processing  System  (IDSIPS)  can  automatically  derive  the  sea- 
surface  slope  statistics  from  sun  glitter  images  through  analysis 
of  the  imaging  geometry.  As  a  late  practical  example,  Fisher  [1986] 
analyzed  four  sun  glitter  images  taken  from  the  space  shuttle 
Challenger  (STS-41G)  to  locate  acoustically  important  oceanographic 
features  in  support  of  hydro-acoustical  sensor  placement. 
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2.2  Spatial  Spectrim  Analysis  of  Naterwave  Scenes 

A  comprehensive  review  would  include  the  %fork  of  Stillwell  and 
Pilon  [1974],  and  Kasevich  et  ai.  [1971,1972].  The  emphasis  of  this 
earlier  work  is  on  the  analysis  of  coherent  optical  processing 
techniques  as  applied  to  photographic  emulsions  of  wave  scenes. 
Kasevich  [1975]  provides  a  general  introduction  to  the  first-order 
theory  subsequent  to  the  development  of  an  geometric-optics- 
approximation  second-order  theory  that  estimates  the  optimum 
viewing  geometry  for  obtaining  reasonable  spectra.  Stillwell  [1969] 
provides  additional  development  of  optical  imaging  theory 
subsequent  to  performing  an  optical  euialysis  to  derive  directional 
wave  spectra.  Only  the  first-order  theory  is  reviewed  here; 
Kasevich 's  second-order  theory  is  outside  the  scope  of  this  study 
since  it  assumes  simplified  approximations  for  both  Fresnel 
reflectivity  and  sky  radiance  distributions.  However,  the  most 
general  results  of  a  second-order  theory  can  be  compared,  with 
caution,  to  the  results  of  this  study;  any  second-order  theory 
loses  definiteness  because  the  spatial  distributions  for  natural 
radiance  are  independent,  i.e.,  no  general  solution  can  be 
specified  [ Stillwell, 1969 ] .  This  is  the  prime  motivation  for  the 
synthesis  and  analysis  of  geometric  effects  through  the  use  of 
analytic  models. 
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Review  of  Kasevich  [1975] 


The  essential  requirement  for  the  determination  of  spectra 
from  wave  images  is  to  have  the  spatial  modulation  of  the  image  be 
proportional  to  the  wave  profile.  Kasevich  uses  the  example  of  a 
transparency  with  film  exposure,  E,  defined  over  its  linear  region 
by 


where 

E(y)-‘fo(y)*f(y) ,  [2.2:2] 

the  mean  exposure  on  film, 

the  exposure  modulation  due  to  scattering  of  radiance 
from  specular  wave-slope  facets,  and 

the  film  gamma.  This  example  is  given  for  a  two- 
dimensional  case. 

If  fo(y)  »  f(y)»  then  Equation  [2.2:1]  can  be  expanded  in  a 
binomial  series  to  yield  the  approximation 


such  that 

fo  (y)  - 
f(y)- 


E{y)-[fAy) 


JL  „ 

nx[l-^ 


f(y) 


2  foiy) 


] 


[2.2:3] 


13 


The  estimation  of  the  slope  spectrvun  from  the  forward  Fourier 
transform  of  the  image  requires  that  f(y)  be  linear  with  respect  to 
the  wave  slope  dz/cty,  where  z  is  the  surface  elevation.  This 
condition  can  only  be  approximately  satisfied  because  of  the 
nonlinearity  of  1)  the  spatial  radiance  distributions  found  in 
nature,  2)  the  Fresnel  reflectivity  variation  with  respect  to 
incidence  angle,  and  3)  the  refraction  of  upwelling  subsurface 
radiance  in  the  direction  of  the  observer. 


Review  of  Stillwell  [1969] 


For  a  small  wave-slope  angle  6,  the  small  angle  approximation 

is: 


p-tan"^  ( 


~  dz 
"  dy 


[2.2:4] 


where  B  is  the  fundamental  parameter  for  extracting  trave-slope 
spectra  from  imagery.  The  law  of  Nalus  defines  the  radiance 
observed  at  azimuth  angle  e  as  a  simple  function  of  Fresnel 
reflectivity  and  incident  radiance: 

Lo(0,  P,«) [2.2:5] 
where 

the  observed  reflected  radiance. 
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the  incident  radiance  to  be  reflected, 

the  Fresnel  reflectivity  (for  any  arbitrary 
polarization) , 

the  zenith  angle  of  observation, 
the  slope  of  the  reflecting  surface  facet, 
the  zenith  angle  of  the  incident  radiance,  and 
the  angle  of  incidence,  such  that 

[2.2:6] 

[2.2:7] 

Figure  2.2:1  illustrates  the  angular  relations  in  two  diaensions. 


FIGURE  2.2:1.  Angular  Relations  for  the  2-D  case. 


1.(1*)- 

/?(«)- 

6- 

P“ 

I*- 

«-6-P 

and 

p-Q-P 
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The  variation  of  observed  radiance  with  respect  to  a  change  of 


surface  slope  dB  at  sone  point  is 


dL(6)  dL(tt)  dp. 
“dp  d|i  dp 


dR(<j) 

du 


do 

dP 


[2.2:8] 


If  fi  is  a  small  angle  and 


then 


~  ♦!?(«)  +L(|i)  *R'(a)  ]  ♦ 

dp 


[2.2:9] 


[2.2:10] 


where  the  prime  (')  denotes  the  first  derivative  with  respect  to 
the  argument.  With  the  assumption  that  the  water  siurface  remains 
analytic,  the  small-angle  linear  approximation  holds  in  the  Fourier 
transform  for  wave-slope  angles  (B)  up  to  at  least  30* 

[Stillwell, 1969]. 

Review  of  Chapaan  and  Irani  [1981] 

The  intent  of  Chapman  and  Irani  was  to  quantify  the  parametric 
dependencies  of  errors  inherent  in  linearly  relating  a  wave-slope 
magnitude  spectrum  to  the  corresponding  radiance  image  magnitude 
spectrum.  Their  approach  was  to  simulate  radiance  images  of  wind- 
driven  water-wave  surfaces  in  three  dimensions.  Their  simulation 
utilized  empirical  models  of  the  water  surface,  the  sky  radiance 
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surface  slope  to  radiance.  They  detemined  that  the  three* 
dieensional  simulation  was  computationally  expensive  so  only  a 
small  number  of  geoMtries  %fere  investigated.  However,  they 
initially  used  a  simpler  t%«o-dimensional  model  of  a  ±15*  sinusoidal 
wave  propagating  in  a  single  direction  along  the  sensor  field  in 
order  to  survey  various  geometries  for  further  synthesis  and 
analysis  in  three  dimensions. 

Figure  2.2:2  illustrates  the  three-dimensional  methodology  for 
image  synthesis  and  error  analysis  used  by  Chai»an  and  Irani. 

Chapman  and  Irani  then  applied  their  three-dimensional  model 
to  six  selected  geometries.  Their  results  indicated  1)  nonlinear¬ 
ities  are  largely  independent  of  mivenunber,  2)  nonlinearities 
remain  relatively  low  (error  <  3  dB)  over  a  relatively  wide  range 
of  azimuth  (100-120*)  aroiind  the  sensor  azimuth,  and  3)  nonline¬ 
arities  increase  within  30-40*  of  the  cross-field  direction. 
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FIGURE  2.2:2.  3-D  methodology  of  Chapman  &  Irani  [1981] 
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Review  of  Korth  [1989] 


North  aeplified  and  extended  the  work  of  Chapaan  and  Irani 
[1981]  based  on  a  proposed  stateaent  of  work: 


"An  exhaustive  evaluation  of  our  radiance  aodulation 
aodel  would  Involve:  1)  siaulating  the  geoaetry  of  a 
surface  area;  2)  transforang  that  surface  geoaetry  into 
an  iaage  using  a  particular  set  of  optical  conditions;  3) 
estiaating  the  slope  spectrua  of  the  siaulated  surface 
and  the  spectrua  of  the  iaage;  4)  repeating  steps  1-3  for 
all  other  possible  iaaging  geoaetries;  and  5)  coaparing 
each  iaage  spectrua  with  the  corresponding  slope  spectrua 
to  obtain  estiaates  of  the  error  arising  froa 
nonlinearity.  This,  in  fact,  is  iapractical  because  there 
are  a  large  nuaber  of  paraaeters  to  be  varied,  and,  for 
any  set  of  these  paraaeters,  the  error  coaputations 
involves  a  large  nuaber  of  operations."  [Chapaan  & 
Irani, 1981] 


The  intent  of  this  study  was  to  develop  and  execute  a  test  design 
«diich  would  provide  a  ainiaal  yet  sufficient  set  of  paraaeters  to 
both  verify  the  liaited  results  of  ChSESuin  &  Irani  and  coaplete  a 
aore  detailed  paraaetric  surface  exploration.  The  aore  robust  test 
design  of  this  study  recognized  that  tdiat  was  coaputation-ally 
iapractical  in  1981  was  quite  doed>le  with  1989  resources. 


Figure  2.2:3  illustrates  the  three-dimensional  methodology  for 
iaage  synthesis  and  error  analysis  used  by  North. 


A  nonlinear  transformation  of  of  wave  slope  to  reflected  and 
refracted  in  both  horizontal  and  vertical  polarizations  was 
effected  under  the  special  conditions  of  Brewster-angle  viewing 
under  clear  skies  at  a  wavelength  of  460  nm.  Seventy  two  syntheses 
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were  varied  with  respect  to  six  distinct  solar  positions,  four 
distinct  wind  directions,  and  three  distinct  wind  velocities. 

The  synthetic  wave  scenes  tiere  analyzed  via  the  forward 
Fourier  transformation  and  their  radiance  magnitude  spectra  were 
compared  with  the  original  slope  magnitude  spectra  of  the  initial 
synthesis  in  order  to  estimate  the  effects  of  the  nonlinear 
radiance  transformation  on  the  recovery  of  the  wave  slope  spectrum 
from  imagery.  Within  the  boundaries  of  the  study,  the  limited 
results  of  Chapman  and  Irani  [1981]  were  more  generally  verified, 
the  existence  of  an  optimal  imaging  geometry  for  slope  spectrum 
estimation  was  indicated,  and  sub-resolution  wave  slopes  were  found 
to  produce  a  significant  effect  on  the  wave  slope  spectra  derived 
from  imagery. 
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3.0 


APPROACH 


The  proposed  research  project  will  have  three  phases:  one 
phase  is  data  synthesis  (modeling  &  prediction);  one  phase  is  data 
collection  (in-field  photography  &  image  processing);  and  one  phase 
is  data  analysis  (of  the  measured  imagery  &  comparison  with  the 
synthetic  model  predictions). 

3.1  Data  Synthesis 

The  3-D  simulation  model  of  North  [1989]  will  be  extended  to 
consider  a  specific  imaging  system:  differentially  polarized  4-lens 
35-mm  panchromatic  photography;  the  photographs  are  digitally 
scanned  to  produce  TIFF  (or  similar  format)  digital  images;  the 
digital  images  are  combined  to  form  the  Stokes'  pareuaeter  images; 
and  these  images  are  transformed  into  their  spatial  spectra  via  the 
forward  2-D  Fast  Fourier  Transform  (FFT).  This  model  will  attempt 
to  simulate  the  actual  imaging  process  that  will  be  executed  during 
the  Data  Collection  and  Data  Analysis  phases. 

A  parametric  exploration  of  the  model  will  then  be  exercised 
in  order  to  predict  the  polarized  imaging  conditions  that  optimize 
spatial  information  content.  The  primary  figure-of -merit  for  making 
this  determination  is  the  amount  of  correlation  between  the  actual 
(modeled)  waterwave  slope  spectra  and  their  reconstruction  from  the 
polarized  radiance  image  spectra.  Other  beneficial  figures-of -merit 
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■ay  be  considered  during  this  phase  of  research. 


One  secondary  objective  (resources  permitting)  is  to  compare 
the  results  of  this  model  with  another  simulation  model,  e.g., 
SENSAT,  using  the  same  input  parameters.  A  correlation  of  results 
between  the  t«K}  models  would  provide  an  added  measure  of  confidence 
in  the  model  predictions  that  will  be  ceunried  into  the  later 
research  phases. 

3.2  Data  Collection 


Differentially  polarized  4-lens  35-mm  panchromatic  photography 
of  various  natural  water  surfaces  will  be  collected  in  the  field 
under  natural  illumination  conditions.  Imaging  geometries,  number 
of  redundant  images,  and  requirements  for  control  and  baseline 
images  will  be  determined  from  the  predictions  of  the  Data 
Synthesis  phase,  above. 

A  Nishika,  Nimslo,  or  similar  4-lens  camera  will  be  modified 
to  produce  the  polarized  images.  One  option  is  to  mount  a 
polarizing  filter  over  each  of  the  four  lenses,  each  filter  rotated 
at  a  different  angle  about  the  optical  axis,  e.g.,  o',  45*,  90*, 
and  135*.  Another  option  is  to  mount  only  three  polarizing  filters 
over  three  lenses,  e.g.,  0*,  60*,  and  120*,  and  use  the  fourth  lens 
to  create  an  unpolarized  calibration  image. 
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The  key  advantage  of  using  a  single-body  4-lens  camera  is  that 
the  synoptic  polarized  images  are  exposed  on  the  same  roll  of  film: 
the  sensitometry  is  equivalent  for  all  the  four-image  sets  on  each 
film  roll.  However,  a  calibration  step  wedge  will  be  imaged  on  each 
roll  of  film  in  order  to  establish  a  Density- versus -log (Exposure) 
curve  for  within-roll  and  between-roll  calibration.  This  exposure 
calibration  is  important  to  the  accurate  measurement  of  signal-to- 
noise  and  contrast  which,  in  turn,  are  key  inputs  to  feature 
recognition  and  reconstruction  algorithms  (e.g.,  surface  slope 
spectrum  reconstruction) . 

The  photographs  will  be  digitized  by  an  image  scanner.  One 
option  is  to  coxamercially  digitize  the  photos  onto  compact  disks, 
e.g.,  via  the  Kodak  Photo  CD”*  process.  A  second  option  is  to 
procure  a  35-mm  slide  scanner,  e.g. ,  the  Nikon  Coolscan*"  LED 
scanner,  and  digitize  the  images  for  direct  transfer  onto  computer 
memory,  hard  disk,  or  other  media.  These  digital  images  will  be 
maintained  as  TIFF,  or  similar  format,  data  files.  The  advantage  of 
the  second  option  over  the  first  option  is  direct  calibration 
control  of  the  digitalization  process,  to  include  the  selection  of 
the  scanning  resolution. 

One  option  for  investigation  is  to  apply  an  additional 
narrowband  spectral  filter  over  the  lenses.  Restriction  to  red- 
bandpass  imaging,  for  example,  should  reduce  the  upwelling 
subsurface  radiance  component  in  the  image;  blue-bandpass  imaging. 
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on  the  other  hand,  should  enhance  it. 


One  secondary  objective  (resources  permitting)  is  to  collect 
polarized  photography  of  waterwave  scenes  containing  Secchi  disks 
suspended  at  various  water  depths  under  varying  conditions  of 
txirbity.  This  secondary  investigation  would  provide  useful 
calibration  data  on  the  upwelling  subsurface  radiance  that  is 
present  in  the  polarized  imagery;  as  an  independent  water 
penetration  study,  it  could  demonstrate  the  utility  of  polarized 
imagery,  compared  with  vinpolarized  Imagery,  in  the  detection  of 
subsurface  features  through  the  differential  reduction  of  unwanted 
polarized  reflected  surface  radiance. 

Another  secondary  objective  (resources  permitting)  is  to 
collect  polarized  photography  of  shadowed  areas  (high  solar  and/or 
high  sensor  zenith  angle)  containing  a  shadowed  target.  This 
independent  investigation  could  demonstrate  the  utility  of 
polarized  imagery,  compared  with  unpolarized  imagery,  in  the 
detection  of  shadowed  features  through  the  differential  selection 
of  polarized  sky  radiance  reflecting  from  the  shadows. 

3.3  Data  Analysis 

The  synoptic  TIFF  images  are  digitally  combined  to  form  the 
Stokes'  parameter  images.  These  images  are  then  transformed  into 
their  spatial  spectra  via  the  forward  2-D  Fast  Fourier  Transform 
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(FFT).  If  required,  nultiple  sequential  realizations  of  the  spectra 
nay  be  conbined  to  provide  statistical  sted>illty.  One  option  for 
investigation  would  first  FFT  the  unconbined  TIFF  images  and  then 
evaluate  combinations  (sum,  difference,  ratio)  of  their  2D  spectra. 

The  majority  of  the  analysis  is  expected  to  involve  the 
evaluation  of  those  polarized  image  spectra  that  demonstrate 
significant  deviations  in  their  distribution  of  spectral  energy 
relative  to  the  unpolarized  image  spectra  for  the  same  imaging 
conditions:  measure  their  ranges  of  spatial  frequency  and  imaging 
azimuth,  measure  their  spectral  contribution  to  the  total  variance 
of  the  scene. 

The  polarized  spectra  of  the  measured  imagery  will  be  compared 
with  their  equivalent  synthetic  spectra  to  determine  the  degree  of 
correlation  between  the  synthetic  and  measured  imagery.  It  is 
anticipated  that  substantial  analysis  may  be  devoted  to  a 
normalization  or  reconciliation  between  the  two  sets  of  results. 
The  outcome  of  a  reconciliation  may  require  modification  of  the 
simulation  model,  re-calibration  of  the  imaging  system,  and  /  or 
discussion  of  any  irreconcilable  differences  and  their  possible 
sources . 

The  expected  study  results  are  the  correlation  sensitivity 
analysis  of  the  simulation  model  data,  reconciliation  with  the 
actual  image  data,  a  comparison  of  analytic  methods,  an 
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Identification  of  methods  that  provide  optimum  results,  any 
unexpected  results,  a  recapitulation  of  the  study  objectives,  and 
finally,  suggestions  for  further  research. 
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4.0 


SCHEDOLB 


MILESTONE 

NOT  EARLIER  THAN 

Siibmit  Proposal 

Spring  Start  94 

Complete  Model 

Summer  Start  94 

Analyze  Model 

Summer  Start  94 

Collect  Data 

Mid-Summer  94 

Analyze  Data 

Mid-Fall  94 

Complete  Thesis 

Mid-Spring  95 

Defend  Thesis 

Spring  End  95 

NOT  IAT8R  THAN 

Spring  End  94 

Suamer  End  94 

Sunmer  End 

Mid-Fall  94 

Mid-Spring  95 

Mid-Summer  95 

Summer  End  95 
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